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Abstract 

Decay rates and time-dependent and direct CP asymmetries in the decays B° — > K + K~ K^^l) 
and KsKsKg^ are studied. Resonant and nonresonant contributions to the three-body decays are 



carefully investigated. Nonresonant effects on 2-body and 3-body matrix elements are constrained 



by QCD counting rules. The predicted branching ratios are consistent with the data within the theo- 
retical and experimental errors, though the theoretical central values are somewhat smaller than the 
experimental ones. Owing to the presence of color-allowed tree amplitudes in B° —* K + K~ Kg^, 
Ph! this penguin-dominated mode may be subject to a potentially significant tree pollution and the 

Qh. deviation of the mixing-induced CP asymmetry from that measured in B —* J/ifsKs, namely, 

A sin 2(3 K + K - Ks(L) = sm.2(3 K + K - Ks{L) — sm.2(3j/^ Ksl can be as large as O(0. 10). In contrast, the 
KgKgKg^ modes appear theoretically very clean in our picture with negligible theoretical errors 
in A sin 2I3k s k s k S ( L) ■ Direct CP asymmetries in K + K~ K S ( L ^ and KsKsKg^ modes are found 
to be very small. 
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I. INTRODUCTION 



Considerable activity in search of possible New Physics beyond the Standard Model has recently 
been devoted to the measurements of time-dependent CP asymmetries in neutral B meson decays 
into final CP eigenstates defined by 

T( j*\ ^ f) ~ T [ m \ - f ] = S, sin(Amt) + Af cos(Ami), (1.1) 
T(B(t)^ f) + T(B(t)^ f) f K J f V ; ' 

where Am is the mass difference of the two neutral B eigenstates, Sf monitors mixing-induced 
CP asymmetry and Af measures direct CP violation (in the BaBar notation, Ct = —Af). The 
time-dependent CP asymmetries in the b — > sqq penguin-induced two-body decays such as B° — ► 
(0, u>, 7T°, rf, fo)Ks measured by BaBar 0, Q] and Belle 0, Q, show some indications of sizable 
deviations from the expectation of the SM where CP asymmetry in all above-mentioned modes 
should be equal to Sj/^ Ks = 0.687 ± 0.032 with a small deviation at most 0(0.1) 0, 0|- Based 
on the framework of QCD factorization 0, the mixing-induced CP violation parameter Sf in the 
seven 2-body modes (<fi, uj, p°, rf, rj, it , fo)Ks has recently been quantitatively studied in 10] and 
It is found that the sign of AS/ = —rjfSf — Sji^k s (rjf being the CP eigenvalue of the 
final state /) at short distances is positive except for the channel p°Ks- After including final-state 
rescattering effects, the central values of AS/ become positive for all the modes under consideration, 
but they tend to be rather small compared to the theoretical uncertainties involved so that it is 
difficult to make reliable statements on the sign at present [loT |. 

Time-dependent CP asymmetries in the b — ► sqq induced three-body decays B° — ► K + K~Ks 
and KsK$Ks have also been measured by B factories BQ,Q,E3,Q,EJJ3 (see Table [J). Three- 
body modes such as these were first discussed by Gershon and Hazumi |l7| |. While KsK$Ks has 
fixed CP-parity, K + K~ K$ is an admixture of CP-even and CP-odd components, rendering its 
CP analysis more complicated. By excluding the major CP-odd contribution from <f>Ks, the 3- 
body K + K~Ks final state is primarily CP-even. A measurement of the CP-even fraction /+ in 
the B° -* K+K-Ks decay yields /+ = 0.89 ± 0.08 ± 0.06 by BaBar 2] and 0.93 ± 0.09 ± 0.05 by 
Belle Q, while the CP-odd fraction in K + R-K L is measured to be /_ = 0.92 ± 0.33±g;^| ± 0.10 
by BaBar ^jj. Hence, while T)f = 1 for the KsK$Ks mode, T]f = 2/+ — 1 for K + K~ K$ and 
rjf = — (2/_ — 1) for K + K~Kl. It is convenient to define an effective sin 2(5 via Sf = —rjf sin2/3 c fj. 
The results of sin2/? e g for K + K~K$ obtained from the measurements of S K + K - Ks and /+ are 
also shown in Table QJ 

In order to see if the current measurements of the deviation of sin2/3 e fj in KKK modes from 
sm.2(3j/ 1 px s signal New Physics in b — > s penguin- induced modes, it is of great importance to 
examine and estimate how much of the deviation of sin2/3 e ff is allowed in the SM. One of the 
major uncertainties in the dynamic calculations lies in the hadronic matrix elements which are 
nonperturbative in nature. One way to circumvent this difficulty is to impose SU(3) flavor symmetry 
Q Kli or isospin and U-spin symmetries j2o| to constrain the relevant hadronic matrix elements. 
While this approach is model independent in the symmetry limit, deviations from that limit can 
only be computed in a model dependent fashion. In addition, it may have some weakness as 
discussed in 
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TABLE I: Mixing-induced CP asymmetries — Sf (top), direct CP violation Af (middle) and branch- 



ing ratios (in units of 10 , bottom) for B — > K + K K$ and KsKsKs decays. For effective sin 2(3 
for K + K~Kq, the third error is due to the uncertainty in the fraction of CP-even contributions to 
the decay rate. Experimental results are taken from iflflQGJEB 



Final State 


BaBar 


Belle 


Average 


K+K-K s a 
(sin 2(3 eG ) K+K - Ks 

K + K-K L h 
(sin 2(3 eS ) K+K - Kh 
KsKsKs 


0.42 ±0.17 ±0.03 
0.55 ± 0.22 ± 0.04 ± 0.11 
0.07±0.28lg;^ 
0.09 ±0.331^1 ±0.10 
0.63±g;|| ± 0.04 


0.52 ±0.16 ±0.03 
0.60 ± 0.18 ± 0.041^2 

0.58 ± 0.36 ± 0.08 


0.47 ±0.12 

„ r „+n 18 

0.07 ±0.30 
0.09 ±0.37 
0.61 ±0.23 


K+K-K s a 
K + K-K L h 
KsKsKs 


-0.10 ±0.14 ±0.04 
-0.54 ± 0.22±g;g| 
0.10 ±0.25 ±0.05 


-0.06 ±0.11 ±0.07 
0.50 ± 0.23 ± 0.06 


-0.08 ±0.10 
-0.54 ±0.24 
0.31 ±0.17 


K+K~K S 
KsKsKs 


11.9 ± 1.0 ±0.8 
6.9±g;| ± 0.6 


14.2 ± 1.7 ±2.0 
4.2^| ±0.8 


12.4 ± 1.2 
6.2 ± 1.2 c 



°with 0(1020)^ excluded. 
6 with (j){l020)K L excluded. 
c with the error enlarged by a factor of S 



1.4. 



We shall apply the factorization approach in this work as it seems to work even in the case of 
three-body B decays 21]. By using factorization and kaon time-like form factors extracted from 
the e + e~ — > KK process, the predicted B° D^ + K~K° rate agrees well with the data |21| . In 
general, three-body B decays are more complicated than the two-body case as they receive resonant 
and nonresonant contributions and involve 3-body matrix elements. Nonresonant charmless three- 
body B decays have been studied extensively |22J, [23J, [2J, [25|, |26J, |27( based on heavy meson chiral 
perturbation theory (HMChPT) 28J, |29J , 30]. However, the predicted decay rates are in general 
unexpectedly large. For example, the branching ratio of the nonresonant decay B~ — > tt + tt^tt^ is 



predicted to be of order 10 5 in |22j and [23j , which is too large compared to the BaBar's preliminary 
result (0.68 ± 0.41) x 10" 6 0. The issue has to do with the applicability of HMChPT. In order 
to apply this approach, two of the final-state pseudoscalars have to be soft. The momentum of the 
soft pseudoscalar should be smaller than the chiral symmetry breaking scale A x ~ 830 MeV. For 
3-body charmless B decays, the available phase space where chiral perturbation theory is applicable 
is only a small fraction of the whole Dalitz plot. Therefore, it is not justified to apply chiral and 
heavy quark symmetries to a certain kinematic region and then generalize it to the region beyond 
its validity. In order to have a reliable prediction for the total rate of direct 3-body decays, one 
should try to utilize chiral symmetry to a minimum. Therefore, we will apply HMChPT only to 
the strong vertex and use the form factors to describe the weak vertex 32]. Moreover, we shall 
introduce a form factor to take care of the off-shell effect. 

As shown in 10], among the aforementioned seven neutral PKs modes, only the ujKs and 
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p Kg modes are expected to have a sizable deviation of the mixing-induced CP asymmetry Sf 
from Sj/^ Ks . More precisely, it is found AS^Ks = °- 12 -ao6 and AS P °K S = -° m -om 1 in the 
absence of final-state interactions Although the tree contribution in these two modes is color 
suppressed, the large cancellation between 04 and a% penguin terms renders the tree pollution 
relatively significant. Unlike the above-mentioned case for two-body decays, the tree contribution 
to the 3-body decay B° — > K + K~ K$ is color- allowed and hence it has the potential for producing 
a large deviation from sin 2(5 measured in B — > J/ ipKs ■ We shall see in this work that it is indeed 
the case. In contrast, the absence of tree pollution in K$K$Ks renders it theoretically very clean 
in our picture. 

The layout of the present paper is as follows. In Sec. II we apply the factorization approach to 
study B° — > K + K~Ks and KsK$Ks decays and discuss resonant and nonresonant contributions 
in Sec. II. Numerical results for decay rates and CP- violating parameters Sf and Af and discussions 
are presented in Sec. III. Sec. IV contains our conclusions. 



II. FORMALISM FOR CHARMLESS 3-BODY B DECAYS 



In the factorization approach, the matrix element of the B — > K KK decay amplitude is given 

by 

(KKK\H eS \B) = ^= E *p(KKK\T p \B), (2.1) 

V ^ p=u,C 

where X p = V p bV ps and 

T p = ai5 pu (ub) V -A ® {su)v_a + a 2 5 pu (sb) V -A &> (uu)v-A + a 3 (sb) V - A &> J2(qq)v-A 

q 

+al^2(qb)v-A ® (sq) V -A + a 5 (sb) V - A <S> E("?<?)v+A 
q <? 

3 

- 2( 4^2(q b )s-P ® (sq)s+P + a 7 (sb) V - A ® E ^ e q(Ql)v+A 

q q 1 

3 3 
- 2a l^2(q~ b )s-P ® -e q {sq) s +p + a g (sb) V - A <g> E i^q^V-A 

q q 
3 

+ a io ^((i b )v-A ® -e q (sq) V -A, (2.2) 
q 1 

with (qq')v±A = 97/^(1 ± 75)?') (Q1 f )s±P = 9(1 ± 75)9' an d a summation over q = u,d,s being 
implied. The matrix element (KKK\j®j'\B) corresponds to (KK\j\B)(K\j'\0), (K\j\~B)(KK\j'\0) 
or (0\j\B)(K KK\j'\0), as appropriate, and a, are the NLO effective Wilson coefficients. In this 
work, we take 

ai « 0.99 ±0.37i, a 2 « 0.19 - 0.1H, a 3 « -0.002 + 0.004i, a 5 ss 0.0054 - 0.005i, 



1 Note that since K + K~K$ is not a pure CP eigenstate, we define Asin2/3 ff = sin2/3 c ff — sin2/3j/^ 
with sin2/3 c ff = —Sf/r]f. In general, the relation A5 1 / = Asin2/3j ff holds for the final state with fixed 
CP parity. 
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o, 4 



-0.03 - 0.02*, a c A 



-0.04 - 0.008i, a v & 



-0.06-0.02*, a, 
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-0.06 - 0.006?, 



a 7 ps 0.54 x 10~ 4 ?', a£ » (4.5 - 0.5*) x 10" 



(4.4-0.3?) x 10~ 4 , 



(2.3) 



0.9 



-0.010 - 0.0002?,, a% « (-58.3 + 86.lt) x 10~ 5 , a c w ra (-60.3 + 88.8?') x 10" 

for typical at the renormalization scale p, = m.b/2 = 2.1 GeV which we are working on. 

Applying Eqs. (|2.1j) . (|2.2I) and the equation of motion, we obtain the £?° — ► K + K~K° decay 
amplitude as 

(K°K+K-\T p \B) = (K+K°\(ub) v „ A \B }(K-\(su) v ^ A \0) [ ai 5 pu + a\ + a{ - {a% + a p s )r x ] 
+{K \(sb) V -A\B )(K + K-\(uu)v-A\0)(a 2 5 pu + a 3 + a 5 + a 7 + a 9 ) 



a 3 + a 5 - + a 9/ 



03 + «4 + 05 - 77(07 + a 9 + a? ) 



-a 
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+{K \(sb)v-A\B°)(K + K-\(dd) v ^ A \0) 

+(K \(sb)v-A\B°)(K + K-\(ss)v-A\0) 

+ (K \sb\B }(K+K-\ss\0}(-2a p 6 + ag) 

+(K + K-K°\(sd) v ^ A \0)(0\(db) v _ A \B ) 

+(K + K-K \s l5 d\0)(0\d l5 b\B°)(-2a p 6 + ag), (2.4) 

with r x = 2m 2 K /(rrib'm s ). In the factorization terms, the KK pair can be produced through a 
transition from the B meson or can be created from vacuum through V and S operators. There 
exist two weak annihilation contributions, where the B meson is annihilated and a final state 
with three kaons is created. Note that the OZI suppressed matrix element (K + K~\(dd)v-A\0) is 
included in the factorization amplitude since it could be enhanced through the long-distance pole 
contributions via the intermediate vector mesons such as p° and u. 

To evaluate the above amplitude, we need to consider the B — > KK, — > KK and — > K KK 
matrix elements, the so-called two-meson transition, two-meson and tree-meson creation matrix 
elements in addition to the usual one-meson transition and creation ones. The two-meson transition 
matrix element (K°K + \(ub)v-A\B°) has the general expression 33] 

{K°(pi)K + (p 2 )\(ub)v-A\B°) = ir(pB-pi-p2)ii, + iu+(p2+Pi)iJi + iv-(p2 

+he flual3 p B {p 2 +pi) a {p2-Pi) 13 - 

This leads to 



■Pi 



(2.5) 



(K~(p 3 )\(su) v . A \0)(K (p 1 )K + (p 2 )\(ub) V -A\B ) 



f_K 

2 



2m|r + (m 2 B - s\ 2 - m 3 l )u + + (s 23 - si 3 -m 2 + m\)oj- 



(2.6) 



where Sj,- = (pi +pj) 2 . A pole model calculation of the B — > K K + transition matrix element 
amounts to considering the strong interaction B° — ► K°B* followed by the weak transition B* — ► 
-RT + and the result is |32| 



(^-(p3)|(^)y-A|0)(K°(pi)K + (p 2 )|(n6)y_ A | J B C 



/x 9^rn B m B 



U s 23 - mi, 



^F{s 23 ,m Bt )F( 



B S K , 



pole 
«23 



m B H m B 

m B 



m B ~ s 23 







(ml)' 



TP B « K ( rn 2\ 

b l (m 3 ). 
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■> . > . (s23 -ml + ml){m 2 B - s 23 - m\) 
rrii + m^ — S13 H 



2m B * 



(2.7) 



where g is a heavy-flavor independent strong coupling which can be extracted from the recent CLEO 
measurement of the D* + decay width, g = 0.59 ± 0.01 ± 0.07 Q, and F B f K are the B s -> if weak 



transition from factors in the standard convention 35j. Since can be far from the mass shell, 
it is necessary to introduce a form factor F(s23,mB*) to take into account the off-shell effect of 
the B* pole. Following |2fJ, it is parameterized as F(s23,ms*) = (A 2 — m 2 B ,)/(A? — S23) with the 
cut-off parameter A chosen to be A = tub* + Aqcd- 

It is worth making a digression for a moment. In principle, one can apply HMChPT twice 
to evaluate the form factors r, uj + and cj_ j^]. However, this will lead to too large decay rates 
in disagreement with experiment [^J. This is because the use of HMChPT is reliable only in 
the kinematic region where K + and K° are soft. Therefore, the available phase space where chiral 
perturbation theory is applicable is very limited. If the soft meson result is assumed to be applicable 
to the whole Dalitz plot, the decay rate will be greatly overestimated. Therefore, we employ the 
pole model to evaluate the aforementioned form factors. We shall apply HMChPT only once to 
the B°K°B* strong vertex and introduce a form factor to take care of the momentum dependence 
of the strong coupling. 

The resonant pole contributions to the form factors r, u± and h can be worked out from 
Eq. I|2.7|) . In principle, there are also nonresonant contributions to these form factors. It turns out 
that the leading nonresonant contribution can be determined as follows. We notice that the same 
B — > KK two-meson transition matrix element also appears in the decay B~ — ► D°K°K~ under 
factorization 21]. The data favors a 1~ configuration in the K°K~ pair j^j. The corresponding 
two-meson transition matrix element is dominated by the u>- term. Following j^] we shall include 
a nonresonant contribution to w_ parametrized as 

U NR = K 2_PB_P1 : (2 . 8) 

and employ the B — > D°K°K~ data and apply isospin symmetry to the B — > KK matrix elements 
to determine the unknown parameter k. The denominator in the above parametrization is inspired 
by the QCD counting rule which gives rise to a l/sf 2 asymptotic behavior, 2 while the numerator 
Pb • P2 = msEx+ is motivated by the observation that K + contains an energetic u quark coming 
from the b — * u transition. 

The matrix elements involving 3-kaon creation are given by 



(K {J (p 1 )K + (p2)K-(p 3 )\(sd) v ^ A \0)(0\(db) v ^ A \B {J ) « 0, (2.9) 



{K\pi)K + {p2)K~{p,)\s lb d\Q)m,bm = 1 - f^L_^_^3 F KKK {m 2 




2 As explained in 21], at least two hard gluon exchanges are needed: one creating the ss pair in K K + , 
the other kicking the spectator to catch up with the energetic s quark to form the K meson. This gives 
rise to a l/s^ 2 asymptotic behavior. 
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where 



m 



ill 



K 



mt 



(2.10) 

m u + m s m s - m d 

characterizes the quark-order parameter (qq) which spontaneously breaks the chiral symmetry. 
Both relations in Eq. (|2.9j) are originally derived in the chiral limit and hence the quark 
masses appearing in Eq. ()2.1U|) are referred to the scale ~ 1 GeV . The first relation reflects 
helicity suppression which is expected to be even more effective for energetic kaons. For the second 
relation, we introduce the form factor F KKK to extrapolate the chiral result to the physical region. 
Following we shall take F KKK (q 2 ) = 1/[1 - (<? 2 /A 2 )] with A x = 0.83 GeV being a chiral 
symmetry breaking scale. 

We now turn to the 2-kaon creation matrix element which can be expressed in terms of time-like 
kaon current form factors as 



(PK + -PK-)„F q K+K - 
? K°K° 

g 



{K + (p K+ )K-{p K -)\q ltl q\Q) 
(K (p K o)K°(p RO )\q 7fl q\0) = (p K o - p K „ ),, Ff " ^ (2.11) 

can be related to the kaon electromagnetic (e.m.) 



The weak vector form factors F^ +K and F^°* y 
form factors F^ K and F^^° for the charged and neutral kaons, respectively. Phenomenologi- 
cally, the e.m. form factors receive resonant and nonresonant contributions and can be expressed 
by 



? K+K~ 



F p + F u + F (j} + F NR , 



F 



K°K° 



-F + F UJ + Fs + F' A 



AT Ji- 



lt follows from Eqs. (J2~TT|) and $ZJM that 



(2.12) 



? K+K~ 



r K+K- 



I'd 
K°K° 



F 



F 



K+K- 



F p + SF^ + -(3F NR 
—F p + 

-3^-^(3*^ + 2^), 



(2.13) 



where use of isospin symmetry has been made. 

The resonant and nonresonant terms in Eq. (|2.12|) can be parametrized as 



F h (s 23 ) 



Ch 



mi 



s 2 3 - im h F h ' 



F (>) 

NR 



/J') J')' 



S23 



'23 



ln|2* 
A 2 



(2.14) 



with A ~ 0.3 GeV. The expression for the nonresonant form factor is motivate d by the asymptotic 
constraint from pQCD, namely, F(t) — ► (l/i)[ln(i/A 2 )] -1 in the large t limit |3j|. The unknown 
parameters Ch,Xi and x\ are fitted from the kaon e.m. data, giving the best fit values (in units of 
GeV 2 for c h ) 

c p( i450) = 7.98 x 10- 3 , 
C w (1650) = -0.116, 



2: 

3c, 



-oj(1420) 



= 0.363, 
-7.64 x 10~ 2 , 



c p(1700) 
c ?i(1680) 



1.71 x 10~ 3 , 
-2.0 x 10~ 2 , 



(2.15) 



and 



xi = -3.26 GeV 2 , x 2 = 5.02 GeV 4 , 



,r-| 



0.47 GeV 2 , 



x 2 



(2.16) 
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Note that the form factors F pui a, in Eqs. ()2.12|) and ()2.13[) include the contributions from the vector 
mesons p(770), p(1450), p(1700), w(782), w(1420), w(1650), 0(1020) and 0(1680). It is interesting 
to note that (i) the fitted values of cy are very close to the vector meson dominance expression 
S v 9vkk for V = p, u), <fr 3£,0], where g v is the e.m. coupling of the vector meson defined by 



(^|jero|0) = gv^v 3 and 9VKK is the V -> KK strong coupling with, -g^x+K- - SpK+K-fV? = 
9uK+K~ /V% — 3.03, and (ii) the vector-meson pole contributions alone yield F^ s +K (0) m 1,-1 
and F^ +K (0) ~ as the charged kaon does not contain the valence d quark. 4 The matrix element 
in the decay amplitude relevant for our purpose then has the expression 

(K ( Pl )\(sb) V - A \B }(K + (p 2 )K-(p 3 )\(qq)v-A\0) = (*i 2 - s 13 )F 1 BK (s 23 )F q K+K - (s 23 ). (2.17) 

We also need to specify the 2-body matrix element {K + K~\ss\0) induced from the scalar density. 
It receives resonant and non-resonant contributions: 

f i—>KK 
^ mf - s 23 - irriiYi 

fs H = ~{ZF NR + 2F' NR ) + v-? r T 

6 s 23 

where the scalar decay constant fi is defined in (z|ss|0) = rriifi, g l ^ KK is the % — > KK strong 
coupling, and the nonresonant terms are related to those in F S K+K through the equation of 
motion. 5 The main scalar meson pole contributions are those that have dominant ss content 
and large coupling to KK. It is found in |4]| that among the fo mesons, only /o(980) and 
/o(1530) have the largest couplings with the KK pair. Note that /o(1530) is a very broad 
state with the width of order 1 GeV |4J|. To proceed with the numerical calculations, we use 



ln(^ 

A 2 



(2.18) 



g f a m)^KK = L5 GeV) 5 /o(i530)-_i^ = 3.i 8 GeV, r /o(980) = 80 MeV, r /o(1530) = 1.160 GeV Q, 
//o(980)(m = rn b /2) ~ 0.39 GeV |42j and // (i 5 30) - // (980)- The sign of the resonantterms is fixed 
by f^ +K (0) = v from a chiral perturbation theory calculation (see, for example, |43j). It should 
be stressed that although the nonresonant contributions to fg and F^ K are related through the 
equation of motion, the resonant ones are different and not related a priori. To apply the equa- 
tion of motion, the form factors should be away from the resonant region. In the large S23-region, 
the nonresonant contribution dominated by the I/S23 term is far away from the resonant one. In 
contrast, the I/S23 term dominates in the low S23-region where resonant contributions cannot be 
ignored. The l/s 23 term in F s is not necessarily conveyed to fs through the equation of motion. 



3 The vector meson e.m. couplings are given by — egm^f^, g pl — [(e„ — ed)/V^]m p f p and g W7 = 
[(e u + ed)/ V^m^f^ where e q is the quark's charge and fy is the vector decay constant. 

4 The sign convention is fixed by using (M(qq' ,p)M(qq' ,p')\qj^q\0) — (M (qq' , p)\qj p q\M (qq' , —p')) — (p — 
p')ti\F q MM \ in the case of a real F q MM . 

5 The use of equations of motion also leads to 

f K+K~ = _ vF K+K- (219) 

Note that the pole contribution to F^ +K should be dropped in the above relation as it applies only to 
nonresonant contributions. 
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Hence, the l/s^ term in Eq. (|2.18l) is undetermined and a new parameter a, which is expected 
to be of similar size as x 2 , is assigned and will be determined later by fitting to the data. The 
corresponding matrix element is now given by 



(K"(p 1 )\sb\B {J )(K + (p 2 )K-(p 3 )\ss\0) 



m 



B 



m K rpBK 
^0 



( S 23)/f + *>23). 



(2.20) 



m b - m s 

Collecting all the relevant matrix elements evaluated above, we are ready to compute the am- 
plitude A(B° -» K S{L) K + K~) = ±A(B° -> K K + R-)/V2. Since under CP-conjugation we 
have K s {pi) -» K s (-pi), K+{p 2 ) -» K~{-p 2 ) and K~{p 3 ) -» K+(-p 3 ), the W -» K S K+K- 
amplitude can be decomposed into CP-odd and CP-even components 

v4[P° -» ^ 5 (pi)K + (p 2 ) J PC-(p 3 )] = A(s 12 , sis, s 23 ) = ^cp- + ^cp+, 

(2.21) 



-4CP± = g [^( S 12> S13, S23) ± ^(S13, Sl 2 , S 2 3)]- 



Correspondingly, we have 

r 



rcp+ + Tcp 
1 



(2vr) 3 32m| 



-^t- [ \A C p±\ 2 dsi 2 dsi 3 
m% J 



1 



/ |A C P±| 2 dsi 2 (is 2 3. (2.22) 



(2vr) 3 32m| 

The vanishing cross terms due to the interference between CP-odd and CP-even components can 
be easily seen from the (anti)symmetric properties of the amplitude and the integration variables 
under the interchange of S12 «-> si 3 - Similar relations hold for the conjugated B° decay rate f . The 
CP-even fraction / + is defined by 



u = 



Tcp+ + Fcp-\ 



r 



(2.23) 



4>Ks excluded. 



Note that results for the K + K~ Kp mode are identical to the K + K~ K$ ones with the CP eigen- 
states interchanged. For example, results for (K + K~ Kl)cp+ are the same as those for 
(K + K~ K$)cp- an d hence /+ in K + K~ K$ corresponds to /_ in K + K~ Kp. 

We next turn to the B° — ► KsKsKs, KgKgKp decays. The decay amplitudes are given by 



A[B° - Ks(pi)K s ( P2 )Ks,l(P3 



1x3/2 

2 



±A[B ^K ( Pl )K (p 2 )K (p 3 )] 
±A[B° -> K\p 2 )K G (p 3 )K°{p l )\ 
+A[B° - K (p 3 )K ( Pl )K (p 2 )}\, 



(2.24) 



with 



K\ Pl )K\p 2 )K\p 3 )\ 



Gp 



^ £ aj (^(pij^^Kdbjv-Ais 

V2 p=uc I 



+(K\p 1 )K\p 3 )\(db) v ^ A \B Q )(K°{p 2 )\{sd) V - A \Q) 



)(K»(p 3 )\(sd) V - A \0) 



(«6 - ^«s) r x) 



+ 



(K U (p 2 )|^|P> U )(^ U (pi)^ U (P3)|^|0) 
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+(K (p 3 )\sb\B )(K ( Pl )K (p 2 )\ss\0)] {-2al + ag) 
^K\ Pl )K (p 2 )K (p 3 )\s l5 d\0)(0\d l5 b\W)(-2al + og) 



+ 



(K»(p 2 )\(sb) V - A \B»)(K»( Pl )K»(p 3 )\(s S ) V - A \0) 



+{K (p3)\(-sb)v-A\B ){K ( Pl )K (p 2 )\(ss) v . A \0) 



a-3 + <H + 05 



(07 + a 9 + aio, 



(2.25) 



where the last term will not contribute to the purely CP-even decay B — > KsK$Ks- Decay rates 
for the KsKgKs and K$K$Kl modes can be obtained from Eq. (|2.22f) with an additional factor 
of 1/3! and 1/2!, respectively, for identical particles in the final state. 

We now consider the CP asymmetries for B° — > K + K~ Kg^, KsKsKs(l) decays. The direct 
CP asymmetry and the mixing induced CP violation are defined by 



Arkk 



Skkk,cp± 
Skkk 



r-r 
r + r 

/ \A\ 2 ds 12 ds 23 - J \A\ 2 ds 12 ds 23 
J \A\ 2 ds 12 ds 23 + / \A\ 2 ds 12 ds 23 ' 

2 J lm{e- 2 ^A CP± A* CP± )ds 12 ds 23 
I \A C p±\ 2 ds 12 ds 23 + f \A C p±\ 2 dsi 2 ds 23 ' 

2jlm{e~ 2lP AA*)ds 12 ds 23 
J \A\ 2 dsi 2 ds 23 + / \A\ 2 dsi 2 ds 23 
/+ Skkk,cp+ + (1 — /+) Skkk,cp-, 



(2.26) 



where A is the decay amplitude of B° — > K + K -Ks(l) or KsKsKs(L)- F° r the K + K K$ mode, 
it is understood that the contribution from </>-?Cs is excluded. It is expected in the SM that 
Srkk,cp+ = sin2/3 eff « sin 2/?, S K kk,cp- ~ - sin 2/? and hence Skx^ « -(2/+ - l)sin2/3. 6 



III. NUMERICAL RESULTS AND DISCUSSIONS 

To proceed with the numerical calculations, we need to specify the input parameters. For the 
CKM matrix elements, we use the Wolfenstein parameters A = 0.825, A = 0.22622, p = 0.207 and 
77 = 0.340, corresponding to (sin 2(3) ckm = 0.724 44]. For B — > K form factors we shall use those 
derived in the covariant light-front quark model 45] with the assigned error to be 0.03, namely, 
F^f(0) = 0.35 ± 0.03. The parameter k in Eq. is determined from the B~ -> D°K°K 



data. From the measured branching ratio B(B~ D°K°K~) = (5.5 ± 1.4 ± 0.8) x 10~ 4 |3J, 
we obtain « = 3.1±f;| GeV where use of of KK = 0.935 and a^ KK {^ a 2 p ) = 0.4 ± 0.2 has been 
made [21]]. For the quark masses and the unitarity angle 7, we shall use m^m^) = 4.2 GeV, 
m s (m b /2) = 80±20 MeV and 7 = (58.6±7)° 44]. The K S K S K S rate sensitive to the parameter a 



6 Writing the CP-conjugated decay amplitude as A = Acp+ + Acp-, we have Acp± = ±Acp± with 
A p -> A*. This leads to S K kk,cp- ~ Skkk,cp+- 
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TABLE II: Branching ratios for B — > K + K Ks, KgKgKs, KsKsKl decays and the fraction of 
CP-even contribution to B° K+K-K s , f+ [see Eq. (JT23J)]. The branching ratio of CP-odd 
K + K~ Ks with 4>Ks excluded is shown in parentheses. Results for (K + K~Kl)cp± are identical 
to those for (K + K~ Ks)cp=f- Theoretical errors correspond to the uncertainties in (i) k, (ii) m s , 
Fq K and a (constrained by the KsKsKs rate), and (hi) 7. 



Final State 


£(10- 6 )theory 




e(10- 6 ) cx pt 


K + K~K S 
(K+K-K S ) CP+ 
(K+K-K S ) C P- 

KsK s K s 
K S K S K L 


r-7 QO+8. 38 + 2. 31 + 0. I 

— 1.08 — 1.59 — 0.10 
c- 4C-+5. 29+1.48+0. 05 
°-1°-0. 65-1. 13-0.06 
1 OQ+3. 08+0. 83+0. 04 
-L-oO-0. 43-0. 46-0.04 
(n /|Q+2.98+0.54+0.03\ 
^•10-0.40-0.22-0.03 ) 

input 

c 7/ | +6.02+2.24+0.02 
°- '^-0.88-1.40-0.03 




12.4 ±1.2 
6.2 ± 1.2 




^•theory 
/+ 




/.expt 
J + 


K+K-Ks 


n nn+0. 06+0. 04+0. 00 
u - yz -0. 16-0.08-0.00 




0.91 ±0.07 




^theory 




jcxpt 


K+K~K L 


n nn+0. 06+0. 04+0. 00 
u - yz -0. 16-0.08-0.00 




0.92 ± 0.37 


TABLE III: Mixing-induced and direct CP asymmetries sin2/3 e g (top) and Aj (in 
respectively, in B° — > K + K~ Ks and KsKsKs decays. Results for (K + K~ Kl)cp± 
to those for [K + K~ Ks)cp^- Experimental results are taken from Table I. 


%, bottom), 
are identical 


Final State 


sin 2/3 eff 




Expt. 


(K + K~K S ) < j > K s excluded 

(K+K-K s )cp+ 

(K + K^Kl)^,K l excluded 

K S K S K S 
K S K S K L 


n 7/IQ+0.080+0.024+0.004 
U. 'iy_o.013~0.011-0.015 

n 77n+ - 113 +°- 040 + - 002 
u -' ' u -0.031-0.023-0.013 
n 7/1 q+0. 080+0.024+0. 004 
U. 'iy_o.013-0.011-0.015 
n 7/1 o+0. 000+0.000+0. 007 
u - ' i°-0. 000-0.000-0. 018 
n 7/IO+0.001+0.000+0.007 
u. 1 _o_ ooi-0 .000-0 .018 




^7+ - 18 
u -°'-0.17 

0.09 ± 0.34 
0.65 ±0.25 




A f (%) 




Expt. 


(K + K~Ks)<f>K s excluded 

(K+K-K s )cp+ 

{K + K^Kl)4,K l excluded 

KsK s K s 
K S K S K L 


n 1 fi +0.95+0.29+0.01 
u - iD -0.11-0.32-0.02 
n nn+0.73+0.16+0.01 
u - uy -0.00-0.27-0.01 
n -ic+0.95+0.29+0.01 
u - iu -0.11-0.32-0.02 
n 7/ i+0.02+0.00+0.05 
u - '^-0.06-0.01-0.06 
n 77+O.I2+O.O8+O.O6 
u - ' '-0.28-0.11-0.07 




-8 ±10 

-54 ± 24 
31 ± 17 



in Eq. (|2.18l) is used to determine a = (— 10.4|t|'g)GeV 4 , where the errors include the uncertainties 
in the KsKsKs decay rate, the strange quark mass and the Fq K form factor. 

Results for the decay rates and CP asymmetries in B° — > K + K~ KgiL), KsKsK S (l) are ex- 
hibited in Table [H] and Table IIII1 respectively. The theoretical errors shown there are from the 
uncertainties in (i) the parameter k which governs the nonresonant contribution to the form factor 
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FIG. 1: The K + K~ mass spectra for B° -> R+R-R s decay from (a) CP-even and (b) CP-odd 
contributions. The insert in (b) is for the <f> region. Results for (R + R~Rl)cp± are identical to 
those for (R + R~Rs)cpt- 



U- [see Eq. (|2.8|) ]. (ii) the strange quark mass m s , the form factor Fq K and a [see Eq. (|2.18|) ] 
constrained from the RsRs^S rate, and (iii) the unitarity angle 7. To compute the CP-even 
fraction /+ and sin2/? e fj for R + R~ Rs, we need to turn off the coefficient in Eq. (|2,13j) . As one 
can see from Table HT1 the predicted rates for B° — > R + R~ Rs(L) decays and the CP-even (odd) 
ratio /+(— ) are in accordance with the data within errors, though the theoretical central values 
on rates are somewhat smaller than the experimental ones. Theoretical errors on the branching 
ratios are dominated by the sizable error in k and the uncertainty in the strange quark mass as the 
penguin term a§r x and the parameter v are very sensitive to m s . Note that the second error in rates 
(including the contribution from the uncertainty in a) are constrained from the RsRs^s rate and 
hence are reduced significantly. For the first error, we note that the larger the value of |k| we have, 
the larger rate on CP-odd R + R~ Rs is obtained, leading to a smaller value of f+(R + R~Rs)- 
Since the central value of our f + {R + R~Rs) agrees well with data, k is preferred to be around its 
central value. 

The R + R~ mass spectra of the B° — ► R + R~R$ decay from CP-even and CP-odd contribu- 
tions are shown in Fig. ^ In the spectra, there are peaks at the threshold and a milder one in the 
large rrix+K- region. For the CP-even part, the threshold enhancement arises from the fo(980)Rs 
and the nonresonant f$ +K contributions [see Eq. (|2.18|) ]. while the peak at large m K + K - comes 
from the nonresonant two-meson transition B° — > R + R$ followed by a current produced R~ . Since 
the nonresonant term [Eq. (|2.8|l ] favors a small m^+Ks region, the spectrum should peak at the 
large m^+K- en d- For the CP-odd spectrum the bump at the large m^+K- en d originates from 
the same two-meson transition term, while the peak on the lower end corresponds to the <pR s con- 
tribution, which is also shown in the insert. The full R + R~ Rs spectrum is basically the sum of the 
CP-even and the CP-odd parts. Note that although we include /o(1530)i?s contribution, its effect 
is not as prominent as one may expect from the R~R + R~ spectrum where a large /x(1500)i^ _ 



contribution is found 



46]. 
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For the mixing-induced CP asymmetry in the K + K~ K$ mode, we compute the effective sin 2(3 
in two different ways: In one way, we calculate S with (pK$ excluded in K + K~ K$ and then 
apply the relation S = —(2/+ — 1) sin2/? e ff and the theoretical value of /+ to obtain sin2/3 e fj. This 
procedure follows closely the BaBar and Belle method of measuring the effective sin 2(3. In the other 
way, we calculate S directly for the CP-even K + K^Kg and identify Skkk,cp+ with sin2/3 e fj. As 
for the KgKgKs mode, there is no such ambiguity as it is a purely CP-even state. As shown in 
Table HTH and Fig. the resulting sin2/? e fj is slightly different in these two different approaches. 

The deviation of the mixing-induced CP asymmetry in B° — > K + K~ Kg and K$K$Ks from 
that measured in B —* J/ipKs (or the fitted CKM's sin 2(3 Q), namely, Asin2/3 e g- = sin2/3 c g — 
s\n2(3j j^Ks (CRM)-, is calculated from Table ITTT1 to be 

Asin2/W-* s = 0.06i°;[| (0.02±g;{g), A sin 2/fr s * s Ar s = 0.06±g;gg (0.02±°;°°). t 3 - 1 ) 

Note that part of the deviation comes from that between the measured sm.2(3jj^ Ks and the fitted 
CKM's sin 2(3. The K + K~ K$ has a potentially sizable A sin 2/3, as this penguin-dominated mode is 
subject to a tree pollution due to the presence of color-allowed tree contributions. For the KgKgKs 
mode, the central value and the error on A sin 2(3 are small. 

It is instructive to see the dependence of sin2/3 e g on the K + K~ invariant mass, mx+K- — 
m 23 = \/S23- For the phase space integration in Eq. (|2.26|) . for a given S23, the upper and 
lower bounds of sx2 are fixed. The invariant mass 77123 is integrated from = 7772 + ^3 to 
77723 = m B — m i- When the variable S23 or 77723 is integrated from to a fixed m^ x (of course, 
Trigg < < ^23)1 the effective sin 2(3 thus obtained is designated as sin 2/3 e fj(77723 ax ). Fig. [21 

shows the plot of sin2/3 e g(777^+ x f ^_) versus 777^+ x R ._ for K + K~ Kg- Since there are two different 
methods for the determination of sin2/3 e fr, the results are depicted in two different curves. It is 
interesting that sin 2(3{rn^ x ) is slightly below svo.2(3ckm at the bulk of the m K + K - region and 
gradually increases and becomes slightly larger than sm.2(3cKM when the phase space is getting 
saturated. The deviation A sin 2(3 K + K - Ks arises mainly from the large m K + K ~ region. 

Direct CP violation is found to be very small in both K + K~ K$ and KgKgKs modes. It is 
interesting to notice that direct CP asymmetry in the CP-even K + K~ K$ mode is only of order 
1CT 3 , but it becomes 0.2 x 1CT 2 in K + K~K S with 4>K S excluded. Since these direct CP asymmetries 
are so small they can be used as approximate null tests of the SM. 

Since direct CP violation in charmless B decays can be significantly affected by final-state 
rescattering we have studied to what extent indications of possibly large deviations of the 
mixing-induced CP violation seen in the penguin-induced two-body decay modes from sin 2(3 deter- 
mined from B — ► J/ipKs can be accounted for by final-state interactions [l0(. It is natural to extend 
the study of final-state rescattering effect on time-dependent CP asymmetries to B — > KKK$ de- 
cays. Final-state interactions in three-body decays are expected to be much more complicated 
than the two-body case. For example, the color allowed tree decay B° — ► Ds + .D(*)~ can rescatter 
into a K+K-Ks final state, where we have D { s * )+ -> K+D*°, £>(*)- -» K S D*- followed by a 
D*°D* -> K~ fusion. These dia grams are too complicated and will not be included in this study. 



7 In passing we note that these diagrams could have the effect of increasing somewhat our predictions for 
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FIG. 2: Mixing-induced CP asymmetry sin2/3 e ff (m^+ x ^-_) (see the text for the definition) versus the 
invariant mass m^+ x ^_ for K + K~ K$ with c/>i£s excluded (solid line) and for CP-even K + K~ K$ 



(dashed line). When m 
saturated and sin2/3 e g(m 
K + K~ Kl mode. 



max 
K+K- 



max 
K+K 



approaches the upper limit m# — mK s , the whole phase space is 
.) is reduced to the usual sin2/? en . This result also applies to the 



Nevertheless, we attempt to incorporate final state rescattering effects in a simple way by includ- 
ing resonance contributions to the corresponding kaon pairs in the final sta te 1471] . We note that 
another attempt in this direction has recently been made by Furman et al. |48j |. They considered 
rescattering of irir and KK pairs in the tttt effective mass range from threshold to 1.1 GeV. While 
their predicted direct CP asymmetry is very small, the parameter S is found to be —0.64 or —0.77, 
depending on the set of penguin amplitudes. However, due to the limitation on phase space, the 
calculated branching ratios of order 1 x 10~ 6 for K + K~ Ks and KsKsKs are only small portions 
of the total experimental rates (see TableCJ and, consequently, the predictions of S may be affected 
when the whole phase space is taken into consideration. 



IV. CONCLUSIONS 

In the present work we have studied the decay rates and time-dependent CP asymmetries in 
the decays B° — > K + K~ Kgijs and KsKsK$(l) within the framework of factorization. Our main 
results are as follows: 

1. Resonant and nonresonant contributions to the hadronic matrix elements are carefully inves- 
tigated. We incorporate final state rescattering effects in a simple way by including resonance 
contributions to the corresponding kaon pairs in the final state. Instead of applying heavy 
meson chiral perturbation theory to the matrix element for B — > KK, which is valid only 



the rates of 2>K final states. Although these contributions carry negligible CP-odd (weak) phases, they 
also contribute to the strong phases and hence will tend to dilute our prediction on AS but not necessarily 
on direct CP asymmetries. 
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for a small kinematic region, we consider the resonant contribution from the B* pole and 
nonresonant contributions constrained by QCD counting rules. 

2. Using the KsK$Ks decay rate as an input, the predicted branching ratio of K + K~ K s ^ 
modes and the CP-even (-odd) fraction of B° — ► K + K~ -K"s(L) are consistent with the data 
within the theoretical and experimental errors, though the theoretical central values on rates 
are somewhat smaller than the experimental ones. 

3. Owing to the presence of color-allowed tree contributions in B — ► K + K~ Kgij^, this 
penguin-dominated mode is subject to a potentially significant tree pollution and the de- 
viation of the mixing- induced CP asymmetry from that measured in B — > J/ipKs, namely, 
Asia.2(3 K + K - Ks = sin 2(5 k +k-k S{L) ~ s ^ n ^Pj/^K s ^ can be as large as O(0. 10). The devi- 
ation A sin 2Bir+ k-Kc< „ arises mainly from the large m K + K - region. 

4. The KsKsKg^ mode appears theoretically very clean in our picture: The uncertainties in 
A sin 2/3 e fj are negligible. 

5. Direct CP asymmetries are very small in both K + K~ Kgi^ and KsKsK$(l) modes. 
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Note added: After the paper was submitted for publication, BaBar has presented a Dalitz 
plot study of B° — > K + K~ Kg decays j^. The BaBar results constrain the tree contribution 
(incorporated via Eq. (2.18) in the present work) in rates and, as a result, a small A sin 2Pk+k-k s 
is preferable. 
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